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TV1 The location problem 4s 


Presenters: Roy Nelson, Wout Oordshoorn and Jaap van der Schaaf 


Rotterdam, in The Netherlands, is one of the largest ports in Europe and is 
an important centre of commerce. It has a population of almost a million. 
Although parts of the town were largely rebuilt after the War, much of 
the traditional style of Dutch housing remains. A major feature of the city 
is the river Nieuwe Maas, which divides the city into two roughly equal 
halves, connected by bridges and tunnels. 


The layout of the city causes practical problems when it comes to the 
location of fire stations. Until the 1960s, the fire stations were run by the 
municipal authority, who divided the city into districts, each served by a 
single fire station somewhere inside it. However, this arrangement 
proved to be unsatisfactory. Several fire stations were located near the 
river, which created a barrier to the movement of fire engines, and so each 
fire station was severely restricted in the area it could conveniently cover. 


In the early 1970s, the city’s services were reorganized. The fire 
department became autonomous, and decided to commission a 
mathematical study from the University of Twente to advise them on 
more satisfactory locations for the fire stations. It was required that there 
should be a response time of at most six minutes to any building in the 
centre of the city. Because the bridges and tunnels crossing the river 
Nieuwe Maas are frequently blocked at rush hours, the parts of the city 
north and south of the river were treated as separate problems. 
Calculations established that the nearest fire station should be within 
2; km of any building in the city. But how did this information help to 
determine the optimum locations for the fire stations? 


An algorithm for solving this type of problem was discovered by Nicos 
Christofides and Peter Viola in the early 1970s. Before outlining their 
method, we first consider a much simpler situation. 


In the following network, the vertices A, B and C represent villages and 
the edges represent roads connecting them. The distances between the 
villages are given in kilometres. 


A single fire station is to be built in such a way that the maximum ' 
distance to any village is as small as possible. Clearly, for practical 
reasons, it must be located on one of the roads or in one of the villages. If it 
must be located in one of the villages, then vertex C must be chosen; from 
this vertex, the maximum distance to any other vertex is 5 km, whereas 
for vertices A and B the maximum distance is 6 km. If, however, it can be 
located anywhere on one of the connecting roads, then the best location is 
at the point 1; km from vertex C along the road to vertex A. The maximum 
distance to any village is then 3; km, which is the least possible value for 
any point on the network. 


Now suppose that, instead of locating a single fire station, we require to 
locate several fire stations on the road network. In some problems, we may 
need to locate the fire stations in such a way that the maximum distance 
from any village to a fire station is less than a specified value. For 
example, how many fire stations must we build so that no village is more 
than 3 km from a fire station, and where on the network do we locate 
them? In other problems, we may be given the number of fire stations that 
are to be built, and we are required to determine where on the network to 
locate them so that the maximum distance from a village to the nearest 
fire station is as small as possible. For example, if we wish to build two 
fire stations so that the maximum distance from a village to the nearest 
fire station is as small as possible, where on the network do we locate 
them? 


These two types of problem are closely related, and are examples of the 
two variants of what is known as the location problem. 


Location problem 


Given a network with n vertices, representing places to be served by a 
number of facilities, the location problem involves one or both of the 
following: 


(a) finding the least number of facilities required so that the 
maximum distance between each place and its nearest facility 
does not exceed a given value d; 


(b) finding the best locations on the network for a given number m of 
facilities so that the maximum distance between each place and 
its nearest facility is a minimum. 


Both variants of the location problem can be solved by the algorithm of 
Christofides and Viola. We illustrate the idea of the algorithm with a 
simple example. 


In the following network, the vertices A-F represent villages and the 
edges represent roads connecting them. The distances between the villages 
are given in kilometres. 

a 


The weights on the edges of the 
network often represent distances, 
and that is how they are interpreted in 
this definition and in the following 
text. 


We want to determine how many fire stations are needed if no village 
should be more than 7 km from the nearest fire station? This is a location 
problem of type (a). We also want to determine the best location on the 
network for two fire stations. This is a location problem of type (b). 


The algorithm makes use of the idea of a region of the network. We define 
a region of a network to be a set of points on the network such that the set 
of vertices of the network that can be reached by travelling a given 
distance 6 from any of the points in the region is always the same. We 
denote a region by the combined labels of the corresponding vertices; for 
example, AB denotes the region from which only vertices A and B can be 
reached by travelling the given distance. A region may comprise just a 
single point, part of an edge or parts of several edges. 


The algorithm works by penetrating further and further away from each 
vertex, along the edges of the network, to construct more and more regions. 
It begins with regions consisting of just the vertices of the network. In the 
case of our example, it begins with the six regions A-F consisting of just 
the six vertices of the network. 


At the first stage of our example, we penetrate 1 km away from each 
vertex, along the edges of the network. This gives the six regions, still 
labelled A-F, shown in the diagram on the left below. From this diagram 
we can draw the bipartite graph shown on the right below, in which the 
vertices on the left represent the regions A—F and the vertices on the right 
represent the villages A—F, and in which a region vertex is joined to a 
village vertex if the village can be reached by travelling at most 1 km 
from some point in the region. 


regions villages 
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6 fire stations needed: 
in A, B, C, D, E, F 


5=1km 


Notice that the bipartite graph and the network diagram enable us to 
deduce that, if no village should be more than 1 km from a fire station, six 
fire stations are needed, one located in each of regions A-F. 


Although the bipartite graph above adds little or nothing to the 
information contained in the network diagram in the case of 1 km 
penetration, we shall see that the bipartite graph provides the key to 
solving the location problem. In general, it is constructed by drawing, on 
the left, one vertex per region of the network and, on the right, one vertex 
per place (vertex) of the network, and then joining a region vertex to a 
place vertex if the place can be reached by travelling at most the current 
distance 6 from some point in the region. 


Let us now penetrate 2 km away from each vertex, as shown on the left 
below. 
regions villages 


5 fire stations needed: 
in-A, BC, DEF 


6=2km 


We now still have the six regions A—F, although they are larger than 
before. But we also have a new region — the point EF, midway between 
vertices E and F. This gives a new bipartite graph, shown on the right 
above, with a new region vertex EF on the left. And, as we can see from the 
bipartite graph and the network diagram, if no village should be more 
than 2 km from a fire station, then only five fire stations are needed — one 
located anywhere in each of regions A—D and one at the point EF. 


Let us now penetrate 3 km away from each vertex, as shown on the left 


below. 
regions villages 


AO---------=252@ A 


3 fire stations needed: 
in AB, CD, EF 


We now still have the seven regions A-F and EF, though they are larger 
than before — in particular EF now consists of a segment of the edge joining 
E and F rather than just a point. But we also have two new regions — the 
points AB and CD, midway between vertices A and B and midway 
between vertices C and D respectively. And, as we can see from the 
bipartite graph, if no village should be more than 3 km from a fire 
station, then only three fire stations are needed — one located anywhere in 
the region EF and one at each of the points AB and CD. 


Let us now penetrate 4 km away from each vertex, as shown on the left 
below. 
D regions villages 


3 fire stations needed: 
in AB, CD, EF 


6=4km 


We now still have the nine regions A-F, EF, AB and CD — all larger than 
before and with AB and CD now segments of edges rather than points. We 
also have two new regions — the points AF and BC, midway between 
vertices A and F, and B and C, respectively. 


The bipartite graph is now more complicated. If we are to use it to 
determine the smallest number of fire stations needed so that no village is 
more than 4 km from the nearest fire station, we need to determine a 
minimum number of region vertices on the left that, between them, are 
connected to all six village vertices on the right. To do this, we first 
simplify the graph by excluding from consideration some of the vertices on 
the left. For example, we exclude vertex A on the left, because a fire 
station located at A can serve only village A, whereas one located in 
either of AB or AF can serve both A and another village. In the same way, 
we exclude all the single-letter region vertices A—F on the left. 


Can we exclude any more vertices? We cannot exclude EF, because the 
village vertex E cannot be reached from any other vertex. Similarly we 
cannot exclude CD. But once we have EF and CD, it only remains to cover 
villages A and B, by means of region vertex AB. 


Thus, if no village should be more than 4 km from a fire station, we still 
need three fire stations — one located somewhere in each of the regions 
AB, CD and EF. We can be more flexible in the location of the fire stations 
than in the 3 km case, since the regions AB and CD are no longer just 
points, but we still need three of them. 


Let us now penetrate 5 km away from each vertex, as shown on the left 


below. 
regions villages 


(ou 


3 fire stations needed: 
in AB, CD, EF 


6=5km 


This time, no new regions are created — we just have larger versions of the 
previous ones. Hence the bipartite graph is the same and, if no village 
should be more than 5 km from a fire station, we still need three fire 
stations — one located somewhere in each of regions AB, CD and EF. 


Let us now penetrate 6 km away from each vertex, as shown on the left 


below. 
regions villages 


3 fire stations needed: 
in AB, CD, EF 

(or AEF, B, CD, 

or AEF, BC, D) 


This time, three new regions are created — the points BE and DE, midway 
between B and E, and D and E, respectively, and the point AEF, 6 km from 
A along the edge AE, from which A, E and F can be reached. 
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In this case, if no village should be more than 6 km from a fire station, we 

still need three fire stations. Even if we choose to locate one at AEF, we still In this case, there are several choices 
need two others to serve B, C and D (for example, in regions BC and D orin Of location for the three fire stations. 
regions B and CD). 


Let us now penetrate 7 km away from each vertex, as shown on the left 


below. 
regions villages 
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2 fire stations needed: 
in AEF, BCD 


Now, we have four new regions — the points CE, ABC, ABF and BCD. We 
can again exclude from consideration all the single-letter vertices A—F on 
the left, and we can also exclude AB, BC, CD and EF, as they are covered 
by the vertices ABC, ABF, BCD and AEF. So, if no village should be more 
than 7 km from a fire station, it is now easy to see from the bipartite 
graph that just two fire stations are needed — one located at the point 
BCD and one somewhere in the region AEF. 


This now enables us to answer both our original questions. If no village Note that, since the regions ABC, ABF 
should be more than 7 km from the nearest fire station, then two fire and BCD are only created as points on 
stations are needed, one at the point BCD and one anywhere in the region Penetrating 7 km from the villages, 
AEF. Also, the best locations for two fire stations are at the point BCD and @"Y distance less than 7 km will 
anywhere in the region AEF, with a maximum distance of 7 km froma “T° Ps ak Ps Aue sis. 
village to its nearest fire station. 


Note that, to solve the ‘how many fire stations should there be if no 
village is to be more than 7 km away from one’ problem, you need only 
perform the last stage above — penetrating 7 km from each vertex and 
drawing the corresponding bipartite graph. 


However, to solve the ‘best location for two fire stations’ problem, you 
need to go through the algorithm stage by stage, as above, until you can 
choose just two region vertices from the bipartite graph to cover all the 
village vertices. In the example, we chose to penetrate in 1 km stages, 
although the distance chosen at each stage will depend on the problem. 


At each stage you should choose a distance such that any new region 
created consists of just a single point. 


The algorithm can be summarized as follows. 


Algorithm of Christofides and Viola 
Given a location problem, START by drawing the network. 
Carry out the following steps for each distance 6. 


STEP 1 Penetrate a distance 5 from each vertex along the edges of the 
network. 


STEP 2 Identify the regions created, and label them with the 
combined labels of all the vertices that are within a distance 
5 of some point in the region. 


Draw the bipartite graph consisting of region vertices on the 
left, representing the regions created in Step 2, and place 
vertices on the right, representing the vertices (places) of the 
network. Join a region vertex to a place vertex if the place is 
within a distance 5 of some point in the region. 


Use the bipartite graph to determine the minimum number 
and best locations of the facilities that can serve all the 
places on the network for the distance 5. 


In the case of location problems of type (a), carry out Steps 1-4 just once, 
with 5 equal to the specified distance d. 


In the case of location problems of type (b), begin with 6 = 0 and 
increase 6 in stages in such a way that each new region created at each 
stage consists of just a single point. STOP when the minimum number of 
facilities determined in Step 4 equals the specified number m. 


The algorithm was used successfully in 1972 to locate the fire stations in 
Rotterdam. It was necessary to use a computer since the calculations were 
too difficult to do by hand. It was found that six fire stations were needed 
north of the river, and five south of the river. As a result, two new fire 
stations were built north of the river, and plans were made to replace 
another at a different location. South of the river, two fire stations were 
closed down and three new ones were built at different sites. 


In the years that have elapsed since that study, there have been many 
changes. For example, a distinction was made between the types of 
building in the different areas; in particular, the response time for a non- 
residential area would generally be larger than that for a residential 
area. Secondly, following the advice of the fire department, many of the 
newer buildings were constructed to a much higher standard, so that fire 
would spread less rapidly through them. This meant that in non- 
residential areas the response time could sometimes be increased to ten or 
twelve minutes, thereby allowing priority to be given to fires in older 
residential areas, where the response time was reduced to five minutes. 


By this time the situation had become very complicated. In order to 
investigate whether any fire station could be closed without jeopardizing 
the adequate fire cover of any district in the city, the Rotterdam Fire 
Service employed a firm of consultants from Apeldoorn. They developed a 
computer model to answer such questions and facilitate the forward 
planning of a modern fire service. By putting the response times for each 
point of the city into a database, they could experiment with the effect of 
moving or closing any number of fire stations. The model was much more 
sophisticated than before, allowing for different response times for up to 
four vehicles. If the resulting plan was not acceptable, the consultants 
could go round the modelling cycle again and again, improving the 
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relocation at each stage. In this way an acceptable plan was produced by 
the planners, the effect of which was known in advance. 


So, back in the 1970s, a relatively straightforward network model 
provided a satisfactory solution to the problem as it was then understood. 
By the early 1990s, the situation had changed in a way that demanded 
considerable refinement of the model. In the future, further developments 
may alter our understanding of the problem yet again, and our 
mathematical model will again need to be changed. But, no matter what 
the solution then becomes, it will almost certainly still be based upon a 
network model. 


TV2_ Tilings at the Alhambra \{.2> 


Presenter: Roy Nelson 


At the foot of the Sierra Nevada mountain range lies the Spanish city of 
Granada. In the fourteenth century, the Moorish kings of southern Spain 
built their magnificent palaces on a plateau overlooking the hillside 
where the Arab artisans had their workshops. The area is famed 
throughout the world, and is known as the Alhambra. 


The first impressions of any visitor to the Alhambra are those of beautiful 
courtyards, such as the Court of the Lions and the Court of the Myrtles. 
Leading from these courtyards are rooms containing elaborate carvings on 
the plasterwork and spectacularly intricate tiling patterns on the walls. 
There are hundreds of such patterns, all seemingly different. But are they 
really different, or can we classify them into types in such a way that we 
can appreciate the similarities between superficially different patterns? 


A simple example of a tiling in the Alhambra is the following pattern of 
square tiles in the Mexuar Court. The tiles are tilted at an angle of 45° and 
are of several different colours, but we can see that it is essentially just a 
square grid. 


A 


Note that the pattern is repetitive — it repeats itself when we move it 
through one or more squares horizontally, vertically or diagonally, and by 
varying amounts in other directions as well. It is the particular 
characteristics of such repetitions that allow us to distinguish various 
types of pattern. And if we rotate the pattern through a quarter of a turn 
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— that is, through 90° — we leave it unchanged; we call this a rotation of 
order 4, since performing it four times takes us back to where we started. 


The next pattern is a mixture of square, hexagonal and octagonal tiles. Its 
outline is a more complicated grid, but again it can be built up by repeating 
the square unit shown. And again, it has a rotation of order 4. 


So, although it looks different from the previous one, it turns out to be 
essentially the same type of pattern. 


Our third example is a plasterwork pattern, built up by repeating a fairly 
complicated square unit. This turns out to be a different type of pattern 
from the previous two. To see this, note that the previous ones look the 
same if we reflect them in a suitably placed mirror. However, this pattern 
always changes if we reflect it, no matter where we place the mirror. If 
we think of the original pattern as going clockwise, then the reflected 
pattern goes anticlockwise. So this pattern is of a different type from the 
others. The earlier ones are ordinary square patterns, but this last one we 
can call a handed square pattern, since it changes under a reflection. 
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There are other ways in which tiling patterns can be similar. Consider the 
following two patterns, found respectively in the Court of the Myrtles and 
in the King’s Hall. 


In the first of these many of the tiles have curved edges, while in the 
second all edges are straight. But although they look very different, they 
both have the property that they are unchanged by a rotation through a 
third of a turn — that is, through 120° — about suitably chosen points. We 
call this a rotation of order 3, since performing it three times takes us back 
to where we started. Since the square patterns we saw earlier have no 
rotations of order 3, and since these patterns have no rotations of order 4, 
they are different types of pattern. 
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However, they do have one similarity. They are both unchanged by a 
rotation through half a turn — that is, through 180° — about suitably 
chosen points. We call this a rotation of order 2. 


And if we rotate the second of these patterns through a sixth of a turn — 
that is, through 60° — about suitably chosen points, we leave the pattern 
unchanged. We call this a rotation of order 6. 


14 


It follows from the above discussion that we can classify patterns 
according to the possible rotations about different points and the possible 
reflections that they permit. It is possible to prove that no pattern can 
have rotations of order k, for any value of k other than 1, 2, 3, 4 or 6, and so 
we can restrict our attention to these values. If we also take reflections 
into account, then it turns out that there are just seventeen different types 
of repeating pattern for tilings, and so any repeating pattern must fall into 
one of these seventeen types. Most, if not all, of these types appear in the 
Alhambra. 


We can extend the above discussion in various different ways. One of these 
is to take account of the different colours of the tiles. For example, the 
following tiling pattern has black and white tiles of the same shape, and 
we can investigate just those rotations that preserve, or interchange, the 
colours of the tiles. In fact, this is what the Dutch artist Maurits Escher 
was interested in when he visited the Alhambra in 1936. Shortly 
afterwards he produced his woodcuts of interlocking birds, such as Day 
and night. 
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Finally, we consider rotations through one fifth of a turn — that is, 
through 72°; these are called rotations of order 5. We have already noted 
that repeating patterns of tiles do not permit rotations of order 5; but if we 
drop the requirement that the patterns repeat themselves at regular 
intervals, then parts of the pattern below do indeed allow such rotations. 
The pattern very nearly repeats itself, but not quite. It was discovered by 
the English mathematician Roger Penrose, and is one of several such non- 
repeating patterns that have been found recently. Such patterns have 
proved to be important in crystallography. 
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TV3 Gas pipeline networks passed 


Presenters: Joe Rooney, Barry Comar, Bill Calderwood, Sam Murdock and 
Luuk Wellens 


In 1965, natural gas was discovered several kilometres beneath the floor 
of the southern North Sea. In the programme we visit the Sole Pit gas 
field to look at the Clipper production platform. This platform became 
operational in 1990 and it exports gas to the mainland. Since that time, 
two satellite platforms have come on stream and a third is currently being 
constructed. 


But where did this gas come from? Why is it under the seabed? To answer 
this we must go back about 250 million years to an era when Britain and 
the North Sea were very different from now. At that time, the familiar 
continents of the world were coalesced into a single giant supercontinent 
known as Pangaea. Northern Europe lay in the arid interior of this land 
mass. The shape of Britain would have been unrecognizable, as would its 
position somewhere close to the Equator. Later, Pangaea began to break up 
and the present-day continents started to drift apart. Britain moved 
north, and the North Sea area was flooded. 


So, over a period of millions of years, Britain and the North Sea area 
passed through several geological stages, ranging from desert to 
submerged marine environment. Layers of mineral and organic sediments 
were deposited until they were buried deep beneath the floor of what 
became the present North Sea. Within these sediments, the great 
pressures and high temperatures (created by the weight of the overlying 
strata) produced both oil and gas from the organic material. 


Just before it was first realized that oil and gas might be found offshore, a 
United Nations convention on the Law of the Sea had been set up. As a 
result, all the countries bordering on the North Sea agreed to divide it up 
into national sectors, and they eventually signed the Continental Shelf 
Act. It wasn’t long before each nation mapped out a grid within their 
sectors and issued exploration licenses to a number of companies to explore 
within particular ‘blocks’. 


The Sole Pit gas field is one of several that have been found in the UK 
sector. It is located almost 100 kilometres off the Lincolnshire coast to the 
North East of the Wash. Although gas was discovered here in 1966, it 
was not the first field to be exploited; some of the others (such as Leman, 
Indefatigable and Sean) were closer to land, or had easier geological 
conditions. The difficult Sole Pit area is only now being developed. Shell 
are currently extracting gas from three reservoirs within the Sole Pit field 
(designated Clipper, Barque and Galleon after the sailing vessels). And 
the name of the game is getting the gas onshore to the customer, quickly, 
efficiently and safely. 
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The terminal at Bacton in Norfolk is where much of the gas from the 
southern North Sea is brought ashore and processed. The complex is 
divided into two separate enclosures; the facility dealing with gas from 
the Sole Pit field is run jointly by Shell and Esso. Once the gas is 
processed, it is sold to the British Gas corporation, who themselves have 
an adjacent facility to distribute the gas throughout the country. The 
control room at the Shell/Esso facility manages the supply of gas from the 
fields run by the two companies. Four main pipelines (called sealines) come 
ashore here from the Sole Pit, Leman, BT and SEAN fields. 


In the programme Barry Comar explains some of the reasons for the 
monitoring and control instrumentation. He describes how Shell starts off 
with a gas daily nomination from British Gas which states how much gas 
must be produced. It is Shell’s task to try to produce that required volume 
over the following twenty-four hour period, at a steady flow rate. Legal 
contracts are involved, so the gas has to be produced within defined limits 
of pressure and temperature and the control room instrumentation at 
Bacton helps the company to achieve this. 


He points out four control panels that relate to the four sealines entering 
the Bacton Terminal, and concentrates on the one that concerns the Sole Pit 
facilities. This monitors the gas produced offshore on the Clipper field 
and sent to Bacton by a pipeline. On one particular panel we see the 
pressure that is in the sealine, and we could observe the flow rates 
through the individual streams, as well as their pressures and 
temperatures. 


It is the pipeline network, the way it is structured and designed, that is 
one of the main features of any gas extraction system. Here at Bacton, 
there is just one 24-inch diameter pipe that comes ashore with all the gas 
from the Sole Pit field. But at the other end of this pipe, there is a rather 
more complicated arrangement. 


Because the gas is trapped in porous sedimentary strata rather than in a 
large hollow reservoir, it flows through the strata at a relatively slow 
rate over large distances. So, the size and shape of the gas fields means 
that gas has to be extracted from several different places. Rather than 
build a completely new pipeline to the mainland for each new well, it is 
sensible to connect the wells into a network — a tree structure. 


But individual offshore installations can themselves reach a larger area 
of the field than you might think, because they use a cluster of wells 
fanning out radially for a considerable horizontal distance, with the gas 
from each well being brought to the surface at a single wellhead platform. 
On the Clipper complex the wellhead platform is where the pipes from 
up to thirty different wells can be combined progressively into common 
pipelines. 


We can view this schematically in two ways. One way is as a ‘star’, with 
all the pipes meeting at a single vertex. We would do that if we wanted 
an overview, with the wellhead platform represented as the single 
vertex. But if we wanted to show more detail, we might represent it as a 
binary tree, with each separate well pipe joining the manifold pipe in 
sequence. 


On this particular wellhead, the individual wells are grouped into two 
clusters — one on each side of the platform. The two groups are brought 
together at a T-junction, and there is a single common pipe crossing to the 
adjacent accommodation platform, where Clipper’s control room is 
located. 


So, although we now have the output from up to thirty wells combined 
into one pipeline, that is not the end of the story. There are currently two 
satellite production platforms some distance away from the main Clipper 
complex. This one is called Galleon, and there is another one called 
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Barque PB. Several more of these satellite platforms are under 
construction, and they all send their gas outputs directly to Clipper. 


The programme shows two pipes running up the leg of the platform. One of 
them comes from the Galleon satellite, which is now on stream. The other 
one is planned to link to a new satellite called Barque PL. Several of these 
satellite platforms are unmanned and are controlled remotely from Bacton 
using telemetry and microwave links, whereas Clipper itself is manned 
and forms a local hub for the others. So our tree is now growing. In addition 
to the wells at Clipper itself, there are links to the wells at Galleon and 
Barque PB, and a link to Barque PL. 


But where does all this gas go to after leaving Clipper? After travelling 
about 70 kilometres, the pipeline arrives on the beach at the Bacton 
terminal. Here the gas is processed, together with gas from three other 
fields (Leman, BT and SEAN) which have separate pipelines to the 
terminal. So, there is a large tree structure of pipelines, bringing the gas 
ashore from an extensive system of offshore production platforms spread 
over a considerable geographical area. This tree has its root vertex at 
Bacton. 


There is another tree structure of pipelines distributing the gas to 
consumers. This has its root vertex at the British Gas complex adjacent to 
Bacton. 


These particular interconnections were chosen and constructed for economic 
and practical reasons. But they are not the only possible tree structures. If 
we have, say, four locations represented by four vertices — the Clipper, 
Barque and Galleon platforms, together with the Bacton terminal — and 
we decide to link them into a tree, in how many ways can we do this? 


The first way is to form a star with Bacton at the hub, and long pipelines 
radiating out to the three platforms. But each of the platforms could have 
been chosen as the hub, giving us three more trees. This gives four different 
stars in all. 


Alternatively, instead of forming a star (where three pipelines meet at 
one vertex), we can form a daisy-chain arrangement with the vertices 
connected serially. There are twelve different ways to do this. Six of them 
have Bacton at the end of the chain, whereas the other six have two 
branches radiating out from Bacton. This gives a total of sixteen different 
ways that we can connect the three platforms and Bacton into a tree 
structure; there are no other ways. 


mae 
PAE 


Barque 


Bacton 


Galleon 


Clipper 


With more platforms, the number of possible trees rises rapidly. If there 
is just one extra platform to incorporate (so that five locations are to be 
interconnected in a tree), then there are 125 different possibilities. With 
six locations there are 1296 different trees. In general, for n locations the 
number of different trees is n”~*, since we are essentially counting the 
number of labelled trees with n vertices. This is a consequence of Cayley’s 
theorem. 


These trees differ from each other in the lengths of pipeline they require, 
and in the geographical routes they take. But how do we decide on the 
geographical routeing of our tree? And how do we connect a new vertex to 
an existing tree? These are some of the factors that Shell engineers have 
to take into consideration when they plan the construction of offshore 
installations. The programme presents a discussion on routeing, with a 
team of three engineers at the Shell Lowestoft depot. 


Bill Calderwood begins the session by focussing on the route that the 
conceptual engineers have suggested for joining the new Barque PL 
platform to the existing facilities. They have proposed a route of about 16 
kilometres from the Barque PL platform to the existing Clipper complex. 
They have not chosen a direct straight-line route, but have selected one 
that swings around in a curve. The team discuss the reasons for that 
choice. They consider a much shorter route of about 9 kilometres from 
Barque PL to Barque PB, rather than to Clipper. 


Sam Murdock observes that the designers avoided a very deep area at a 
location called Coal Pit, which the shorter direct route would have to 
cross. The longer route raises no particular problems along its track and it 
coincides with the existing pipeline corridor on the approach to Clipper. 
Also, there are no wrecks along the seabed. However, it is known that 
there may be some sand waves up to about 5 metres high on the seabed in 
an area lying about 3 kilometres from the Barque PL platform. There seem 
to be no serious problems with the longer route. So the choice of route 
cannot be dictated just by the length of pipe. There may be other 
considerations — technical, legal, geographical and economic. 


Bill Calderwood then raises the issue of having to install extra metering 
facilities on the platform to which the new pipeline is joined. There may 
not be the physical capacity to do this, and it might necessitate signing 
separate legal agreements. Luuk Wellens joins the discussion with a 
comment that the existing metering arrangements rule out the building of a 
new metering platform. 


Sam Murdock re-inforces this view by stressing that unmanned platforms 
are to be preferred, because it is then possible to carry out maintenance at 
fewer intervals, rather than having to install additional costly facilities 
on Barque PB for manned operation. Bill Calderwood raises the potential 
problem of the new pipeline having to cross an existing one owned and run 
by another company. Sam Murdock acknowledges this, but is not unduly 
concerned. Negotiating agreements with other companies or other 
developers to cross their pipelines should not pose a problem, but Shell 
would have to look at it. The company would also have to look at the 
areas where Shell must have other third party agreements. However, in 
the case of Barque PL to Clipper, this pipeline stays within the 
boundaries of Shell’s own licence areas. 


So there are lots of complications to be considered, as well as distance. 
Licence areas and pipeline crossings are among the legal constraints, 
seabed terrain is a geographical constraint, gas composition and pressure 
provide technical constraints, and the need to meter the flows down 
shared pipelines can provide an economic constraint. Taking account of 
these constraints may even suggest a structure with cycles. This might 
arise when an older tree structure can take only some of the output from a 
new platform; an additional pipeline would be needed for the remainder 
of the output. 
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Even after a route has been established, the technical complications have 
to be considered. Some of these are dealt with, once the gas reaches the 
terminal at Bacton. For instance, the gas is unlikely to be pure methane, 
but it will contain varying amounts of other gases and liquids such as other 
hydrocarbons, carbon dioxide, water and small quantities of helium. 
These need to be reduced or removed. 


One way this is done is by using what is called a ‘slug-catcher’. This is an 
underground pipe arrangement that changes the flow of the gas, and 
allows ‘slugs’ of liquid impurity to drain off at the bottom. There is one of 
these slug-catchers for each of the four sealines at Bacton. Once the liquid 
has been removed, other impurities such as water vapour can be dealt 
with by reducing the temperature of the gas. This is done by using a 
chiller unit, cooled by liquid propane, so that the impurities can be drawn 
off in the cold separator tanks. 


Another impurity removed here is glycol. This is an antifreeze, 
deliberately added to the gas at the platforms to prevent the formation of 
crystalline hydrates in the pipelines. The glycol can be recycled, and it is 
returned to the platforms in small separate pipelines. Technically, these 
return pipelines introduce cycles into our tree structure, but since they carry 
no gas we can ignore them. 


The one remaining process to be carried out on the gas is filtering. This 
removes any solid particles, so that the final product meets its contractual 
level of quality. The gas can now be analysed and metered before it is fed 
across the road to the British Gas distribution terminal. 


The programme shows how the gas flows through these processes from the 
mimic diagram in Bacton control room. Sealine 4 is the one from Clipper, 
and after passing its slug-catcher, the gas moves through a group of valves 
to the pipe labelled gas stream number 4. From here, it passes through the 
chillers and dust filters to the final metering area. 


At this level of detail, the pipes no longer form a tree structure. There are 
alternative routes through the processing plant in order to maximize 
flexibility. The valves can be set to close off any cycles in the network so 
that the final gas export is indeed the root vertex of a tree. 


Unfortunately, the conditions are not always right for this purification to 
work effectively. In some circumstances, impurities condense within the 
offshore pipeline before reaching Bacton. These deposits constrict the gas 
flow and would eventually block the pipe. 


To remove any blockages, a sphere — known, colloquially as a Pig — is 
sent down the pipeline from Clipper to Bacton. It is pushed along the 
inside of the pipe by the gas pressure, and clears out any deposits. When it 
arrives at Bacton, it is collected in the appropriately-named ‘sphere 
receiver’. 


Of course the pressure of the gas has to be high enough to push the spheres 
along the pipeline; in fact, the gas pressure at the platform is another 
important consideration in the design of the system. 


Initially, the gas leaves the under-seabed reservoir at a high enough 
pressure to be transported directly along the pipelines to the shore. But as 
more and more gas is extracted, the reservoir pressure drops, so compressors 
have to be used to ensure an adequate flow along the pipeline. These are 
just large turbines, typically powered by aircraft engines of up to 30,000 
horsepower. 


Periodically, of course, the compressors need to be taken off-line for 
maintenance. This has to take place without interrupting the gas flow, 
and again this means that there must be alternative routes for the gas. 
This creates yet more cycles when we look at the network in detail. 
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As well as compressing the gas, its flow has to be measured. This 
information is not just needed for monitoring purposes: when gas from 
different platforms is sent down the same pipe, the volumes from each 
platform can be calculated from the flow rates and used to determine the 
composition of the mixture. And if the platforms are operated by different 
companies, then the volumes are also used to determine the correct 
payments for the gas. This also means that some of the gas may leave the 
main pipeline and travel along a parallel route through the transducers, 
leading to yet more cycles in a detailed view of the network. 


Unlike some of its satellites, Clipper is a manned platform with its own 
small control room. Any problems with the operation can be monitored 
here using an array of alarms, and remedial action can be taken where 
necessary. But overall control of the network is carried out remotely. The 
radio masts are not just for verbal communication: they relay telemetry 
signals to and from the terminal at Bacton. And it is from here that valves 
in the pipeline network can be opened and closed to alter the gas flow as 
needed. | 


So, with chillers, slug catchers, compressors and meters, we do not have a 
tree when we look closely at the detailed gas flow. However, when 
looking at the general flow of gas we can aggregate these detailed 
processes into single vertices or edges and treat the overall gas collection 
and delivery system as a tree structure. And since the root vertex is where 
the production companies sell their gas to the customer, you could say that 
money really does grow on trees. 


TV 4 Kinematics 


Presenter: Joe Rooney 


Flying a jet aircraft takes a great deal of skill, and pilot-training is an 
expensive and potentially risky business. So nowadays, instead of using 
real aircraft, many airlines train their pilots on flight simulators. These 
offer a much safer alternative and allow trainee pilots to make what in a 
real situation might be dangerous mistakes without killing themselves or 
others in the process. They can also practise their responses to unusual or 
emergency scenarios and explore the operational limits of their plane. 
Then, having mastered one type of plane, they can use the same basic 
simulator (suitably re-programmed) for further training on a different 
aircraft. 


But how does a simulator reproduce the motion of a real aeroplane in 
flight? Any flight simulator must generate, coordinate and feed back 
forces and visual impressions, so that pilots can experience convincing 
sensations of flight. The simulator must respond to movements of the 
controls, so as to achieve a high degree of realism in motion. This way, 
the pilot’s senses are fooled. Intuitively, it must look and feel just like a 
real flight deck. So it is not surprising that airlines need highly 
specialized equipment. This often includes a standard kinematic sub- 
system, as in the simulator at Hughes Rediffusion shown in the 
programme. This is because all aircraft must be capable of performing the 
same basic types of motion — to a greater or lesser extent. 


To see why this is, let us run through a typical take-off sequence for a real 
commercial airliner. The plane accelerates down the runway until it 
reaches its flying speed. The nose comes up in the so-called ‘rotate’ 
manoeuvre to increase aerodynamic lift on the wings. And the aircraft 
leaves the ground in a tilted-back climbing attitude. It then banks to the 
left or right in order to turn onto its eventual heading, and continues to 
climb quite steeply. Eventually it reaches its allotted cruising height, 
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drops the nose, and aligns the wings into straight and level flight. 
Throughout the flight the plane generally maintains this status, except 
for occasional mid-course corrections to its heading or altitude. Near its 
destination, it begins the landing sequence. The nose drops and the plane 
starts to descend. Further banking manoeuvres then align it with the 
centre of the runway as it makes its final approach to touchdown. Note 
that, although the nose is pitched down for most of the descent, the 
sequence ends with a nose-up attitude to gain maximum lift at the 
relatively slow landing speed. Once on the runway, the engines reverse 
thrust and the aircraft decelerates to a safe taxiing speed. 


If we now look at a military fighter plane, we see essentially the same 
basic take-off and landing sequences, although higher speeds and 
accelerations are achieved. The big difference here is that the aircraft 
can perform more complicated aerial manoeuvres such as vertical climbs 
and dives, 360-degree rolls, loops, and spins. If we look at a Harrier jump- 
jet or a helicopter we see even more unusual aerial ‘choreography’, such as 
spinning around a vertical axis, rapid sideways movements, and even 
flying backwards. 


The motion of each of the three types of aircraft (the jetliner, the fighter 
and the helicopter) involves a continuous change of position and 
orientation in space. So, how do we describe the location of an object in 
three-dimensional space? 


To locate an object in space, we need a set of numbers, or coordinates. For 
the position of the object, three coordinates are needed to give the position 
of, say, the object’s centre of mass, and usually these three numbers 
measure the position along three directions mutually at right angles (such 
as the x, y and z directions, or latitude, longitude, and height above the 
ground). For a moving object such as an aircraft we often refer the 
movement to coordinate directions associated directly with the object’s 
attitude (such as forwards, sideways and upwards). The technical names 
for straight-line movement (also called translation) in these coordinate 
directions are 


surge, for translation in the forwards (or backwards) direction, 
sway, for translation in the sideways (left or right) direction, 
heave, for translation in the upwards (or downwards) direction. 


But the position of the object is only half the story. We need three more 
numbers to describe its orientation, and usually these three numbers 
measure angular rotations about axes (straight lines) in the coordinate 
directions. The technical names for rotational movement about coordinate 
axes attached to the object are 


roll, for rotation about a fore-and-aft axis (along the object), 
pitch, for rotation about a transverse axis (across the object), 
yaw, for rotation about a vertical axis. 


So we need a total of six coordinates (three for position and three for 
orientation) to describe what we call the pose of the object. We say that 
the object has six freedoms. If we want a kinematic system, such as a 
flight-simulator, to be able to move an object into a given pose, then it 
must provide these six freedoms in its motion. In other words, six 
independent inputs must drive and control the system. 


One way to achieve a given pose is to use a robot manipulator arm. The 
diagram illustrates a PUMA robot (Programmable Universal Manipulator 
for Assembly), which can grasp an object in one pose and transfer it to 
another pose along a specified path. If the robot were used in a factory 
environment, the object might be a welding torch for spot-welding car 
components, a paint spray-gun for painting body panels, or an electronic 
control unit being assembled into, say, a washing machine. But if the arm 
were much bigger, it could carry a larger object such as a mock-up of an 
aircraft flight-deck — it would be a potential flight simulator. 
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The PUMA robot is a serial arrangement of links and joints. In fact, there 
are seven rigid links in the system, connected in series by six joints. Each 
joint provides one rotational degree of freedom, and this must be driven 
and controlled independently. The largest joint rotates the whole arm, 
relative to the base. The arm combines the six joint freedoms to give a 
total of six freedoms at the robot’s hand (the gripper). This means that 
the hand can move in the six standard ways — surge, sway, heave, roll, 
pitch and yaw. But in general, several joints must operate simultaneously 
to produce the correct coordinated motion of the hand. For instance, three 
joints must be activated to give a simple straight-line translation of the 
hand. The activation of a single joint does not produce such a simple 
change in the pose of the hand, where only one pose coordinate alters, and 
vice versa. Simple joint movements give complicated movements of the 
hand, and simple movements of the hand require complicated 
combinations of joint movements. 


So why are robot arms not used to simulate flight? Because the robot is a 
serial arrangement of links and joints, any slight error in the positional 
control of a joint angle would be amplified as we move along the arm, and 
the errors in successive joints would accumulate. The greatest 
amplification of errors occurs for the first joint at the base, and the effect 
decreases the closer we get to the hand. The dynamic effects when the arm 
is moving make matters even worse. So the final pose of the hand would be 
difficult to control dynamically. Try writing with your arm and fingers at 
full stretch and you'll appreciate the problem. And notice that, for many 
delicate tasks or for precise control, you often need to use two hands (and 
arms) cooperating together in a non-serial arrangement. 


In contrast to the robot, a typical commercial flight simulator has a non- 
serial arrangement of links and joints. The configuration is based on a so- 
called Stewart platform, first proposed by Stewart in the early 1960s. 
Instead of a hand supported and controlled by the movements of an arm, 
the flight simulator has a platform supported and controlled by the 
movements of six legs. The movements needed to simulate a real flight are 
derived by computer from the pilot’s controls, and they can be very 
complicated indeed. So the programme considers a more basic set-up 
consisting of a stripped-down skeleton version of a commercial flight 
simulator. 


The skeleton rig has the same arrangement and type of legs as a 
commercial flight simulator, but the cockpit is replaced by an equivalent 
weight of steel girders — about eighteen tons’ worth. The legs supporting 
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the flight simulator form a parallel arrangement of links and joints. The 
system consists of fourteen rigid links, connected in parallel by six 
hydraulic rams (the legs) and twelve further joints at the ends of the legs. 


Each ram provides two freedoms — one is a rotation and the other is a 
translation — but only the translational freedom is driven and controlled 
independently. The six rams combine their single translational freedoms 
in parallel, so as to provide the required total of six freedoms at the 
cockpit platform. 


Just like the robot hand, the platform can be moved in six independent 
ways, and the skeleton rig has facilities to select each of these. It can 
surge, sway, heave, roll, pitch and yaw. But again, like the robot arm, 
several joints must be activated simultaneously in order to produce the 
correct coordinated motion of the platform. The six freedoms of the 
platform’s legs must be orchestrated carefully, so as to achieve each 
required change of pose. The activation of a single ram does not produce a 
simple change in the pose of the platform (where only one pose coordinate 
alters), and vice versa. 


Unlike the robot, the six-legged platform is a good kinematic design for a 
flight simulator. Because the platform is a parallel arrangement of links 
and joints, any errors in the control of joint variables are not amplified, nor 
do they accumulate successively; they may even cancel out at the 
platform. So the final pose of the hand is easier to control dynamically. 
But it does have its drawbacks. The parallel arrangement is far more 
restricted in its range of movements than the serial design. 


We can represent the structure of a kinematic system using a graph. The 
graphs for the robot and simulator platform are as follows. 


P26, 


These are interchange graphs; each rigid link of the system is represented 
by a vertex, and each joint is represented by an edge. For the simulator, 
the cockpit platform is represented by a vertex of degree 6 as indicated, 
and the six joints connecting the platform to its six legs are represented by 
six edges, as are the joints connecting the legs to the base on the ground, 
which is itself represented by another vertex of degree 6. We can see from 
the graphs that the platform should be more stable than the robot. The 
multi-cycle graph of the platform has five more paths for distributing 
and balancing the mechanical forces and controls than does the tree graph 
of the robot. You could think of the robot as one-legged (a monopod) like a 
pogo stick, whereas the simulator is six-legged (a hexapod) like an insect. 


The structure of a kinematic system can be modelled by an interchange 
graph, but its actual motion rests on the design of its joints. The PUMA 
robot has six joints, all of the same type, and each joint permits only a 
rotational movement. The simulator platform has eighteen joints, but this 
time there are two types. Each ram permits a combination of rotation and 
translation, where the rotation usually has quite a small range. The other 
joints are called pin-in-sphere joints; these permit a combination of two 
different rotations about non-parallel axes. These joints are quite 
complicated in their design, but their function can be understood more 
simply in terms of some primitive types of joint. 


There are six basic types of joint that allow controlled movement between 
two links in contact at surfaces sliding over each other. They are known as 
the Reuleaux pairs. 


¢ A revolute pair allows a rotation about one axis. This is only a single 
freedom, and so the joint has imposed five constraints on the relative 
motion of the two links. The contacting surface can be any surface of 
revolution. The programme uses the simplest — a cone — generated 
by a straight line revolving around the axis. 


e¢ A prismatic pair allows a translation along a direction. This joint 
also imposes five constraints. The contacting surface can be any 
extruded surface. The programme uses the simplest — a triangular 
prism — generated by three straight lines translating along the axis. 


e A screw pair allows a rotation about an axis, combined with a 
translation along that same axis. This is still just a single freedom, so 
again the joint imposes five constraints. The contacting surface is like 
a screw thread and this is used in the programme, although the 
simplest such surface — a helicoid — is generated by a straight line 
at right angles to the axis, rotating about and translating along the 
axis. 


¢ Acylindric pair allows a rotation about an axis and an independent 
translation along that same axis. This gives two freedoms, so the 
joint imposes only four constraints. The contacting surface is a 
cylinder. 


e A planar pair allows independent translations in two different 
directions and also a translation about an axis perpendicular to these 
two directions. This gives three freedoms in all, so the joint imposes 
only three constraints. The contacting surface is a plane. 


e — Aspherical pair allows independent rotations about three mutually 
perpendicular axes. This joint also imposes three constraints. The 
contacting surface is a sphere. 


The six joints of the PUMA robot function as revolute pairs, even though 
they have a more complicated internal structure with gears, ball- 
bearings, and so on. So, in kinematic terms, the motion of the hand is 
produced by six revolutes. On the simulator platform, the six hydraulic 
rams function as cylindric pairs, and each of the twelve pin-in-sphere 
joints appears to function as a combination of a revolute pair and a 
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spherical pair. But before we look at this in detail, let us see how the 
motion of a system arises from the separate motions of its links and joints. 


The programme shows a simple example to illustrate the procedure. It’s 
an analogue of the three-dimensional PUMA robot arm, but it is two- 
dimensional; all of its movements are constrained within a plane. Each of 
its joints is a revolute, and each joint axis is perpendicular to the plane of 
movement. It consists of four links, connected in a serial arrangement by 
three revolute joints. 


The hand has three freedoms. It can move horizontally and vertically, 
and it can also rotate within the plane. So it can adopt an arbitrary two- 
dimensional pose. We say that its mobility is 3. : 


To calculate the mobility, we start at the first link. That’s the base, so it 
is stationary. 


If the second link were free to move in the plane, it would have three 
freedoms on its own. However, it is joined to the first link and the joint 


constrains its movement — it can no longer move horizontally or vertically . 


and can only rotate about the joint axis. So the joint removes two freedoms: 
there is only one left. We say that the joint imposes two constraints. 


The third link adds three more freedoms — and the second joint removes 
two of them. So this moving link has accumulated a second freedom 
relative to the base link. 


The same happens for the fourth and final link and the third and final 
joint. This gives three freedoms for the hand at the end of the chain. The 
hand therefore has a mobility of 3. 


In the calculation of the mobility, we’ve taken the number of links, 
subtracted 1 — that’s the number of moving links — and multiplied by 3, 
because each moving link would add three freedoms if it weren’t joined to 
any other link. Then we’ve subtracted the number of revolute pairs, 
multiplied by 2, because each revolute pair removes two freedoms. The 
result is the total mobility of the system. If M is the mobility, n is the 
number of links, and j is the number of revolute pairs, the formula is: 

M = 3(n - 1) - 2). 
Another two-dimensional example shown in the programme is an 
analogue of the flight simulator. This has prismatic pairs as well as 
revolute pairs, and the translational movement of the prismatics is 
parallel to the plane. You can see that the platform has three freedoms 
and can adopt a range of different two-dimensional poses. It can move 
horizontally and vertically, and it can rotate within the plane, so its 
mobility is 3. But how is this mobility calculated? 
We proceed as before. In this system there are six links, but only. five of 
them move since the base is fixed in position. So, there are 6 — 1 = 5 moving 
links, each contributing three freedoms. 
The four revolute pairs remove two freedoms each, and the two prismatic 
pairs also remove two freedoms each, since, just like the revolute pair, 
each prismatic pair constrains the joined links to have only one relative 
freedom (translation, in this case), instead of the three freedoms they 
would have if they weren’t joined. 


So the mobility is 
3n(o=t)=4x2=208 
which comes to 3, as expected. 


In general, for a two-dimensional system containing revolutes and 
prismatics, if M is the mobility, n is the number of links, jg is the number 
of revolute pairs, and j, is the number of prismatic pairs, the formula for 
the mobility is: 

M = 3(n — 1) — 2jp — jp. 
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In a three-dimensional system, the motion of the links is no longer 
parallel to any given plane. For example, revolute pairs do not generally 
have their axes parallel, and prismatic pairs do not usually have their 
axes parallel to a single plane. And the four other primitive types of joint 
may be present also; for example, there are cylindric pairs on the 
simulator platform. In general, the axes of these joints are neither 
parallel, nor perpendicular, nor intersecting. They are said to be skew. So 
the motion is much more complicated, and the mobility calculation 
involves several more terms. 


In three-dimensional space, an unconstrained link has six freedoms (three 
for its position and three for its orientation). But a revolute pair joining 
two links allows only a single relative freedom — a rotation about an axis 
— and so it imposes five constraints. So does a prismatic pair — it allows 
only a translation along one direction. And so does a screw pair — it 
allows a combined rotation about, and translation along, an axis, which is 
still just a single freedom. But a cylindric pair is different. It allows two 
freedoms — a rotation about, and an independent translation along, an 
axis — and so it imposes only four constraints. A planar pair allows three 
freedoms — translations in two independent directions, and a rotation 
about an axis perpendicular to both of these directions — and so it imposes 
only three constraints. So does a spherical pair, which allows rotations 
about three mutually perpendicular axes. 


We illustrate this by calculating the mobility of the PUMA robot. There 
are seven links, one of which is stationary, so we start with a total of 6 x 
(7 — 1) = 36 freedoms. There are six revolute pairs, so we subtract 5 x 6 = 30 
constraints, leaving a mobility of 6, as expected. If we control each joint, 
then we can choose an arbitrary three-dimensional pose for the hand. 


In general, consider a three-dimensional system containing all six types of 
primitive joint. If M is the mobility, n is the number of links, j, is the 
number of revolute pairs, jp is the number of prismatic pairs, j,, is the 
number of screw pairs (H stands for helix), j. is the number of cylindric 
pairs, j, is the number of planar pairs (E stands for Euclidean plane), and 
jg is the number of spherical pairs, then the formula for the mobility is: 


M = 6(1t- 1) — 5jp — Sip — Sis ~ Sic Sis — Sis 


Now, a serial robot is suitable for supporting only comparatively light 
loads. We might try supporting a heavier load with a configuration 
which appears to be similar to the arrangement on the flight simulator — 
there is a base, a platform to support the load, and six legs each in two 
parts, giving fourteen links altogether. In this case, the two links in each 
leg are joined by a prismatic pair, so that the length of the leg can be 
varied; this gives six prismatic pairs. The ends of the legs are attached by 
spherical pairs, so that each leg can rotate in three-dimensional space; 
this gives twelve spherical pairs. If we now perform the mobility 
calculation with these numbers, we get: 


M=6x(14-1)-5x6-3x12=12. 


Surprisingly, the mobility formula gives a mobility of 12, whereas we 
expect a mobility of 6, since the system should be capable of controlling 
the three-dimensional pose of the platform. We seem to have six extra 
freedoms. Where are they? What has happened is that each leg has a 
spherical pair at each end, so that rotation is allowed at both ends. This 
means that the leg is free to rotate as a single unit about its own 
(longitudinal) axis; it is called a passive freedom, because it does not affect 
the overall configuration of the system. But it does affect the mobility 
calculation, since the six passive degrees of freedom (one for each leg) in 
the system account for the missing mobility. Passive degrees of freedom 
are undesirable, because they are not normally under direct control. 
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In a real flight simulator, the dynamics of the platform must be totally 
under control, and so passive freedoms have to be avoided. This means 
that spherical pairs cannot be used at both ends of the legs. So what can 
we use instead? We could use a spherical pair at just one end of each leg 
with some other type of joint at the other end, but in practice it is better 
not to use any spherical pairs. 


We have already seen that the real simulator uses pin-in-sphere joints at 
the ends of each leg, so that neither end can rotate about the 
(longitudinal) axis of the leg. But this would remove too much freedom (by 
imposing too many extra constraints) if the legs remained prismatic. The 
non-passive rotation that does need to remain in each leg comes from a 
joint in the centre of the leg — a cylindric pair, rather than a prismatic 
pair — and this allows rotation about, as well as translation along, the 
axis of the leg. This means that each hydraulic ram has two freedoms 
instead of one, and each end has a joint with two freedoms instead of 
three. This eliminates the passive freedom in the leg. 


But what type of joint is this pin-in-sphere joint? How does it relate to 
the six types of primitive pair? And how do we obtain the correct 
mobility value using the mobility formula? In fact, the pin-in-sphere 
joint allows the same motion as a particular type of universal joint, known 
as a Hooke’s joint, named after Robert Hooke who developed it in the 17th 
century from a much older device. This type of universal joint is made of 
two revolute pairs with their axes perpendicular to each other, and so it 
allows two independent rotations. The pin-in-sphere joint can move in the 
same way. The outer sphere can revolve about the pin axis, giving one 
rotation, and the slot in the outer sphere can slide along the pin, giving a 
second rotation about an axis perpendicular to the pin axis. 


So let us analyse a platform using these universal (Hooke’s) joints to 
attach the legs. Each universal joint counts as two revolute pairs, together 
with an extra link — the cross-piece inside the joint. So there are twenty- 
six links altogether — two from each of the six legs, one from each of the 
twelve universal joints, one for the base and one for the platform. There 
are twenty-four revolute pairs (two from each universal joint) and six 
cylindric pairs (one from each leg). Carrying out the calculation, we 
obtain the mobility formula 


M=6x (26-1)-5 x 24-4x6=6, 
giving a mobility of 6, as expected. This parallel arrangement of links and 


joints gives us the six freedoms we want, so that pilots no longer need 
planes for realistic flight training. 
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TV5 The four-colour theorem 


On 13 October 1852, Augustus De Morgan, Professor of Mathematics at 
University College, London, wrote a letter to Sir William Rowan 
Hamilton, the Astronomer Royal of Ireland. 


My dear Hamilton, 

A student of mine asked me today to give him a reason for a fact which I did not 
know was a fact — and do not yet. He says that if a figure be anyhow divided 
and the compartments differently coloured so that figures with any portion of 
common boundary are differently coloured — four colours may be wanted, but 
no more... the following is the case in which four are needed. 


pa de | 


A, B, Cand D are 
names of colours 


My pupil says he guessed it in colouring a map of England. The more I think of 
it, the more evident it seems ... 


This problem later came to be known as the map-colouring problem, or as 
Guthrie’s problem, after the student (Francis Guthrie) who first posed it: 
is it true that any map can be coloured with just four colours so that no two 
neighbouring countries have the same colour? 


Note that countries meeting at a point are not considered to be 
neighbouring countries; for example, country A is a neighbour of both B and 
D, but it is not a neighbour of C — so it is permissible for both A and C to be 
coloured the same, and for both B and D to be coloured the same. 


aa, 


Note also that it is immaterial whether we include the outside region in 
our colouring, since it can be regarded as an extra (ring-shaped) country, as 
shown below. We shall usually omit the outside region. 


blue 


The fact the we need at least four colours, and not three or two, can be seen 
from the map in De Morgan’s letter in which all four countries meet each 
other, and more complicated maps such as the following one. 


In order to consider the problem mathematically, we must specify exactly 
what we mean by a map. We consider a map to be a plane drawing of a 
connected planar graph, drawn so that its faces correspond to the countries 
(apart from the infinite face which is the outside region mentioned 
above). Since the problem is to colour the faces on the two sides of each 
edge differently, we assume that the graph has no bridges and in 
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not allowed 


no bridges no vertices of degree 2 


particular, no vertices of degree 1. We also assume that the graph has no 

vertices of degree 2, since such vertices do not affect the colouring of the 

faces; more generally, we exclude all graphs with cutsets containing two Asin Graphs 3, we express these ideas 
edges. So the vertices of a map have degree 3 or more, and we will usually _ by saying that a map has no cutsets 
omit our conventional black spot vertex label. with less than three edges. 


Definition 


A map is a plane drawing of a connected planar graph containing no 
cutsets with 1 or 2 edges. 


It was not until the late 1870s that any progress was made in the solution. 
In June 1878, seven years after De Morgan’s death and 26 years after the 
problem was first posed, the Cambridge mathematician Arthur Cayley 
asked the following question at a meeting of the London Mathematical 
Society: 

Query: has a solution been given of the question as to whether four colours are 
sufficient to colour all maps? 


Shortly afterwards, Cayley drew attention to a difficulty that lies in the 
way of providing a simple proof: 

In any tolerably simple case, it can be seen that four colours are sufficient, but I 
have not succeeded in obtaining a general proof: and it is worthwhile to explain 
wherein the difficulty consists. Supposing a system of n areas coloured 
according to the theorem with four colours only, if we add an (n+1)th area, it 
by no means follows that we can without altering the original colouring colour 
this with one of the four colours. For instance, if the original colouring be such 
that the four colours all present themselves in the exterior boundary of the n 
areas, and if the new area be an area enclosing the n areas, then there is not any 
one of the four colours available for the new area. 


However, Cayley showed that, in trying to prove that four colours are 
sufficient for all maps, it is sufficient to prove the result for those maps in 
which exactly three countries meet at each intersection point, such as the 
following map. 


In 1879, Alfred Kempe, a London barrister and keen amateur 
mathematician, presented a proof that four colours are always sufficient; 
his paper appeared in the newly-founded American Journal of 
Mathematics. 


Kempe first observed that every map must contain somewhere within it a 
country with at most five neighbours — this can be proved using Euler’s 
formula. Thus every map must contain at least one of the following: 


¢ acountry with just two boundaries — a digon; 

¢ acountry with just three boundaries — a triangle; 
* acountry with just four boundaries — a square. 

¢ oracountry with five boundaries — a pentagon. 
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The exact shapes of the countries do 


not matter: only the number of 
z boundaries is important. 


digon triangle square pentagon 


Since every map must contain one or more of these, we call this set of 

countries an unavoidable set: it cannot be avoided. Similarly, it can be 

deduced from Euler’s formula that every map must contain at least one of . 
the following configurations of countries, so these also form an 

unavoidable set. 


digon triangle square two pentagons pentagon and hexagon 


In proving the four-colour theorem, Kempe assumed that the statement 
that four colours are sufficient is incorrect, and supposed that some maps 
need five or more colours. He then showed that such an assumption leads 
to a contradiction. We now describe Kempe’s argument. 


If there are maps that need five colours or more, then among all such maps 
there must be at least one with the smallest number of countries. Let us 
consider such a map in detail — this map needs five or more colours, but 
any map with fewer countries needs at most four colours. 


We now consider the four configurations in turn. 


Digon If the map contains a digon, we remove the digon by shrinking it 
to a point. The resulting map has one fewer country and so, by our 
assumption, we can colour its countries with at most four colours. What 
happens when we reinstate the digon? Since it has only two neighbours, 
which can be coloured with at most two of the four colours, there is a spare 
colour to colour the digon. This gives a four-colouring of the original map, 
contrary to our assumption, and so yields the required contradiction. 


remove digon . colour map replace digon 


red red . 


Triangle If the map contains a triangle, we shrink it to a point. The 
resulting map has one fewer country and so, by our assumption, we can 
colour its countries with at most four colours. What happens when we 
reinstate the triangle? Since it has only three neighbours, which can be 
coloured with at most three of the four colours, there is a spare colour to 
colour the triangle. This gives a four-colouring of the original map, 
contrary to our assumption, and again yields the required contradiction. 


remove triangle colour map replace triangle 


So that deals with maps containing a digon or a triangle. We say that a 
digon or triangle is reducible because if we reduce either to a point and 
four-colour the remaining map, then we can extend the colouring to a four- 
colouring if the whole map, with the digon or triangle included. 


Je 


Square If the map contains a square, we shrink it to a point, as before. 
The resulting map has one fewer country and so, by our assumption, we can 
colour the countries of this new map with at most four colours. What 
happens when we reinstate the square? Here a difficulty arises, since the 
square has four neighbours which may use all four colours, so there may be 
no spare colour available to colour the square. 


colour map 


At this stage Kempe applied an ingenious argument. He reasoned that we 
can always arrange to colour the four neighbours with just three colours. To 
do this, we choose two of them that are not adjacent, such as the red and 
green neighbours. Note that the red neighbour may interconnect with 
several other red and green countries, and similarly for the green 
neighbour. 


If the two red-green parts are separate from each other, then we can 
interchange the reds and greens in, say, the upper part without upsetting 
the colouring of the other part. 


Now only three colours (green, blue and yellow) are used to colour the four 
neighbouring countries, and so the square can be coloured with the fourth 
colour, red. We can therefore four-colour the whole map, and obtain the 
required contradiction. 


However, if the upper and lower red-green parts are not separate from 
each other, but link up in a chain as follows, then interchanging the reds 
and greens in the chain does not help. 


replace square 


a 


gS 
. 


Berane 
oe se 


lue 
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But now, because of the solid red-green chain, the blue-yellow parts on the 
right and left of the square must be separated from each other, so we can 
interchange blues and yellows in the right-hand part without affecting 
the blues and yellows on the left. As a result, only three colours (red, green 
and yellow) are now used to colour the four neighbouring countries, and so 
the square can be coloured with the fourth colour, blue. So, again we can 
four-colour the whole map, and obtain the required contradiction. 


So the square is reducible — when we remove it, four-colour the remaining 
map, and then reinstate the square, we can always extend the colouring to 
the square. 


Pentagon Kempe applied the same process to the pentagon. In this case, 
the difficulty arises when the pentagon has five neighbours that use all 
four colours. To avoid this, he changed two of the colours around the 
pentagon at the same time. 


By the argument above, we may assume that there is a red-green chain 
from the lower red neighbour to the upper-right green one, and that there 
is a red-yellow chain from the lower red neighbour to the upper-left 
yellow one. 


Thus we can interchange the blues and yellows on the right and the blues 
and greens on the left, so that the only colours next to the pentagon are red, 
yellow and green. The pentagon can therefore be coloured blue. So again 
we can four-colour the whole map, and obtain the required contradiction. 


Since every map contains a digon, triangle, square and/or pentagon, and 
since we have proved the theorem in each case, this completes the proof 
of the four-colour theorem — or so Kempe thought! 


Kempe’s proof was quickly accepted, and became part of mathematical 
folk-lore. For example, Lewis Carroll formulated the problem as a game 
between two players: 

A is to draw a fictitious map divided into countries. B is to colour it using as 
few colours as possible... A’s object ts to force B to use as many colours as 
possible. How many can he force B to use? 


And a Headmaster set the map-colouring problem as a challenge to his 
school: 

it is found, on trial, that four colours are always sufficient, whatever the 
shape or number of the countries or areas may be. Required, a good proof of 
this. Why four? .... Solutions to be sent to the Head Master before Dec. 1st... No 
solution may exceed one page, 30 lines of MS., and one page of diagrams. 


Even the Bishop of London got hooked. He thought he had proved it 
while allowing his mind to wander during a meeting. But all he had 
proved was that there do not exist five neighbouring regions in the plane. 
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This does not prove the result — if five neighbouring regions exist, then 
the four colour theorem is false; but proving that five neighbouring regions 
do not exist does not establish that the four colour theorem is true. 


So the problem seemed well and truly settled. But shortly afterwards, in 
1890, Percy Heawood published a paper Map colour theorem that pointed 
out a flaw in Kempe’s proof of the map-colouring theorem. 


The present article does not profess to give a proof of this original theorem; in 
fact, its aims are so far rather destructive than constructive, for it will be shown 
that there is a defect in the now apparently recognised proof. 


Heawood produced an example that demonstrated where Kempe’s 
argument was incorrect. Recall that in trying to prove that the pentagon is 
reducible, Kempe had carried out two interchanges of colour at the same 
time. Either interchange on its own is fine, but the following coloured map 
presented by Heawood shows that in general it is not possible to perform 
both interchanges simultaneously — the object is to perform the 
interchanges of colour in such a way as to colour the pentagon in the 
middle. 


The blue and yellow countries adjacent to the pentagon are connected by a 
blue-yellow chain, which separates the reds and greens above the 
pentagon from those below it. 


So we can interchange the reds and greens in the upper part without 
upsetting those in the lower part, and we obtain the following colouring. 


ee 


Alternatively, we could have done another interchange. The blue and 
green countries are connected by a blue-green chain, which separates the 
reds and yellows above the pentagon from those below it. 


So we can interchange the reds and yellows in the lower part without 
upsetting those in the upper part, and we obtain the following colouring. 


Either colour interchange is permissible on its own, but what happens if 
we do both at the same time? We get two red countries coming together, as 
follows. 
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originally 
green 


originally 
yellow 


So, in general, we cannot perform two colour interchanges at the same 
time. Thus Kempe’s proof is incorrect. 


This error that Heawood discovered in 1890 highlighted a major 
difficulty. However, two of Kempe’s ideas proved to be crucial to further 
progress: those of unavoidable set and reducible configuration. Let us remind 
ourselves of these. 


We started with an unavoidable set of four configurations — the digon, 
triangle, square and pentagon. We saw that three of these are reducible, 
in that any colouring of the rest of the map can be extended to include this 
configuration as well. But we were not able to cope with the pentagon. So 
we need to replace it by other configurations of countries to see whether 
these are reducible — if not, we try another unavoidable set, and so on. For 
example, if we take the pentagon and replace it by two adjacent pentagons 
and a pentagon joined to a hexagon, then we get another unavoidable set, 
as we saw earlier. By replacing these by yet more configurations, we get 
larger and larger unavoidable sets containing more and more complicated 
configurations. We need to show that we can eventually obtain an 
unavoidable set of reducible configurations. But can such a set be constructed? 


The answer came 124 years after Augustus De Morgan’s letter. The Times of 
23 July 1976 contains the following story. 


Two American mathematicians have just announced that they have solved the 
proposition that has been puzzling mankind for more than 100 years. It ts 
simply that four is the maximum number of colours needed on any map so that 
no two contiguous countries are coloured the same... Their proof published 
today runs to 100 pages of summary, 100 pages of detail and a further 700 pages 
of backup work. It involved 1,000 hours of computer time, 10,000 diagrams 
and the computer printout stands four feet high on the floor. 


The two mathematicians were Professors Kenneth Appel and Wolfgang 
Haken of the University of Illinois. In the television programme and in 


the following extract from the New Scientist of 21 October 1976, they 


describe how they went about constructing an unavoidable set of reducible 
configurations. Since then, thousands of unavoidable sets of reducible 
configurations have been discovered, so that there are now thousands of 
proofs of the four colour theorem. 


Enter the Fast Digital Computer* 


Even though a great deal of progress was made in the study of reduci- 
bility, the goal of proving the four colour theorem by demonstrating the 
existence of an unavoidable set of reducible configurations seemed 
extremely far off. The critical problem was this: No one had any 
reasonable intuition of a set of configurations which was unavoidable and 
seemed to contain configurations which were likely to be reducible. In 


Recall that unavoidable means that 


every map must contain at least one 


of them and reducible means that 


whichever one it is, then the proof can 


be completed. 
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particular, if there were such a set it was not clear that it was small 
enough so that its members could be tested for reducibility. 


With the advent of large, fast digital computers, a new tool was made 
available to workers in the field. Heesch formalised the ideas of Kempe 
and he and his students used computers to show a great many 
configurations reducible. Heesch strongly believed that an unavoidable 
set of configurations could be found and used a method which has since 
come to be known as the principle of discharging (by analogy with the 
idea of moving charges in an electrical network) to find an unavoidable set 
of reducible configurations for maps with certain restrictions. 


Wolfgang Haken, in the late 1960s, noticed that Heesch’s discharging 
arguments could be greatly improved and simplified. He argued that since 
the study of reducibility had proceeded much farther than that of 
unavoidable sets, one should spend a great deal more effort on the study of 
unavoidable sets. From the work on reducibility, especially the many 
configurations studied by Heesch, it became evident that configurations 
with certain easily checkable properties were rather likely to be 
reducible while those without these properties were unlikely to be 
reducible. 


In 1972, Haken and I began to search for unavoidable sets of configurations 
which were likely to be reducible. We did this by examining various 
discharging procedures to see what sort of sets of configurations would be 
generated. A drawback to this approach had been that the study of a 
single discharging procedure would take several months and would not 
help greatly in simplifying the study of a second procedure. We overcame 
this difficulty by creating a very large sophisticated computer program 
which could, by minor variations in input and certain parameters, be used 
to study a great many discharging procedures. This has the advantage 
that while the program took a long while to perfect, once it was 
available it was possible to study the results of a procedure in a few hours. 


By late 1974, we became convinced that unavoidable sets of likely-to-be- 
reducible configurations consisting of a few thousand configurations each 
could be found. It also appeared that the amount of computer time required 
to check the reducibility of the configurations in such sets would be large 
but not prohibitive. Unfortunately, the problem remained that if a single 
configuration in such a set were irreducible the set would not serve the 
intended purpose of proving the four colour theorem. 


It appeared, however, that a method could be developed to modify such 
sets to replace unwanted configurations. Haken, and John Koch, and I then 
combined efforts on developing a collection of computer programs to test 
configurations for reducibility. Although such programs had been written 
before (by Heesch, S. Gill, Allaire and Stewart, and others), because of 
technical differences in purpose it was thought advisable to have a new 
set. 


By January 1976 it appeared that the study of discharging procedures had 
reached a point in which a serious attempt on the four colour theorem 
could be made. The final attempt used a rather flexible technique in 
which one configuration at a time was generated for the potential 
unavoidable set. When it was generated an immediate attempt was made 
to show it reducible. If this could not be done with reasonable effort (up to 
30 minutes on an IBM 370-168 computer) it was discarded and the 
procedure was modified to avoid its use. Previous work had given us 
confidence that this method would converge to an unavoidable set of 
reducible configurations. In June 1976, after analysis of 10,000 
configurations, over 2000 of which were tested for reducibility, and after 
using over 1000 hours of computer time on various computers owned by the 
University of Illinois, an unavoidable set of under 2000 reducible 
configurations was produced. 


(* Reprinted, with permission, from New Scientist, 21 October 1976.) 
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TV6 Compact discs ..6 


Presenters: Roy Nelson, Stan Baggen, Jack van Lint 


When you listen to a compact disc, it is obvious that the sound quality is 
far superior to that of early recording media. This is partly due to the use 
of advanced technology, particularly the use of a laser to scan the disc. 
But there is more to it than that — the secret lies within the compact disc 
itself. 


Each compact disc has a spiral track, like a gramophone record. But, 
unlike the record, the track is scanned outwards from the centre, rather 
than inwards, and it is scanned at constant speed (just over a metre per 
second), so that the disc slows down as it plays, from eight revolutions per 
second down to four. 


Most importantly, the sound is recorded on the disc in a digital format. 
Any sound, whether music or speech, can be represented as an analogue 
signal. With stereo, there are two such waveforms, for the left and right 
channels. To convert each analogue signal into digital form, it is sampled 
at regular intervals, and the value is measured. This value is then 
represented as a binary number — a string of sixteen bits, each 0 or 1. These 
bits are processed eight at a time. Each group is called an information 
symbol, and the symbols are encoded before they are recorded onto the 
compact disc. This combination of laser scanning and digital recording was 
developed in the Netherlands by Philips in the early 1980s. 


In the programme a large-scale model of a compact disc is used to 
demonstrate some of the key features of the system. Light is reflected off 
the under surface of the disc as it rotates. Sometimes it is bright and 
sometimes dark, but it is in the transition from light to dark and back 
again that we find the coded information. 


Along the spiral track are depressions (called pits) in the surface of the 
disc, while the parts of the track level with the surface are called lands. 
The pits have different lengths but the same depth, just less than a 
quarter of a wavelength of the laser light. So the light reflected back 
from the bottom of a pit travels about half a wavelength further, and is 
therefore nearly 180° out of phase with the light from the adjacent 
surface. Destructive interference makes the pits look darker than the 
lands. Each land/pit or pit/land transition corresponds to a 1-bit, with 
0-bits in between at regular intervals. 


To make a recording one needs a master disc, created using a high-powered 
argon laser. The beam is carefully focused onto the blank disc and 
modulated with the digital signal. Once the master disc is completed, the 
process of creating sub-masters and final production discs is then purely 
mechanical. 


There is also much technology involved in reading a disc. The laser spot on 
the surface of the disc has a diameter of about one micron. This has to 
follow the spiral track accurately, even if the turntable is slightly off- 
centre. To show how this can be done, the programme contains an 
experiment using a server system with three spots, one central spot that 
reads the information, and two satellite spots on either side of the track. 
If the track moves away due to an eccentricity, one satellite spot sees less 
of the track, while the other sees more. This imbalance in reflected light 
then causes the server system to correct the position of the central spot 
that reads the information. In the demonstration, the servo system 
controls a mirror that deflects the spot from side to side, so that it can 
follow horizontal variations in the track position. But there is also 
vertical motion — up to a millimetre — which needs to be followed. 


Ch 


With the help of the experiment, we can see how the optimum distance 
between the lens and the disc is maintained. The light is first focused onto 
the disc by the focus servo, and on its return path passes through a prism 
and is focused onto two small spots on these detectors. The positions of 
these two spots tell us the optimum distance between the lens and the disc, 
as follows. If the distance between the lens and the disc changes, those 
two spots on the diodes move inward or outward, depending on the 
distance between the lens and the disc. The servo system then responds by 
moving the lens back to the optimum position with respect to the disc. 
Now the lens is mounted onto a coil that moves in a permanent magnet, just 
like a loudspeaker system. By making use of error signals, as shown 
earlier in the experiment, we can optimally control the distance of the 
lens with respect to the disc. 


So, sophisticated control systems ensure that the laser spot follows the 
track as closely as possible. But so far there is nothing to suggest any 
improvement in the quality of the recorded signal, and this turns out to be 
essential for the success of the system. 


In fact, the compact disc system would not be feasible without an efficient 
error correction system, since cheap mass replication of gigabits of 
information is impossible without creating any errors. Already a freshly 
pressed disc may contain hundreds of thousands of errors — and handling 
the disc may create even more, through fingerprints or dirt, for example. 


So, with digital representation of sound, error correction can be used to 
improve quality. In the CD recording process, the signal is first converted 
into digital form as described above. It is then passed through an error 
correction encoder, and finally it is sent to a special modulation encoder 
which averages out the number of 1s and 0s recorded on the disc. The 
reverse happens when the disc is played back. The signal first passes 
through a modulation decoder to remove the averaging effect. Then it goes 
to an error correction decoder. Finally, it is converted back to analogue 
form to be sent to the amplifiers and loudspeakers. 


But how does error correction work? We have seen that music is recorded 
digitally as a long string of 1s and Os. Now, a error may occur in reading 
the information, so that a 1 is interpreted as a 0, or a 0 may be interpreted 
as al. 


For example, consider the short piece of ‘music’ 1011. Let us send seven bits 
instead of four, and let us write the four bits into the parts labelled 1, 2, 3 
and 4 in the diagram shown in the margin. The remaining three bits are 
determined by making the number of 1s in each circle even; for example, in 
the third circle (at the bottom) there are two 1s and one 0, so the 
remaining bit is 0. This gives the seven-bit message 1011010. 


Now suppose that an error occurs in position 3, so that the 1 in position 3 is 
read as a 0. The receiver sees that something is wrong in circles 1 and 2, 
because in each circle the number of 1s has become odd. Since circle 3 is 
correct, the error must have occurred in the overlap between circles 1 and 2, 
but not in circle 3. Now that we know where the error is, we can change the 
erroneous symbol into the correct one. This code corrects 1 error, and we 
have used just seven symbols instead of the original four. 


The error-correcting codes employed on compact discs use a special type of 
arithmetic, called modular arithmetic. We illustrate this by considering 
the numbers from 0 to 12, which are added or multiplied modulo 13 — that 
is, we add the numbers as usual, and then take the remainder on dividing 
by 13. For example, 8 + 11 is normally 19, but if we take the remainder on 
dividing by 13, we get the answer 6. Similarly, 6 x 8 is normally 48, but if 
we take the remainder on dividing by 13, we get the answer 9. The 
following are the addition and multiplication tables for this type of 
arithmetic. 
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Note that we can also divide in this arithmetic. If we look down the 
column for 4, say, then each number from 0 to 12 occurs just once. For 
example, 9x 4=10, andso10+4=9. 


We now come to the Reed-Solomon code, which is the type of code used on 
the compact disc. For simplicity, we shall not use long strings of 1s and 0s, 
but assume that the music is recorded using the integers from 0 to 12, with 
arithmetic carried out modulo 13. 


Let us assume that the music is given as a string apa, . . . 439, and as before 
we add some extra symbols. The original symbols are cg, cy, .. . , C19, and we 
assume that we have two further symbols cj, and cj, determined by 
requiring that 

Cott, +... + Cy =O, 


and that if we multiply each number by its location and add all of these 
numbers, the sum is again 0 — that is, 


Cy +2cp +33 +... +12c47 = 0. 


These two equations determine the final two symbols c;; and cp. 


Let us now assume that an error occurs in the position with number i, so 
that instead of c;, we receive r; =c; + e, where e is the magnitude of the 
error, and the sum cy +... + C2 is now e, instead of 0. In the second equation, 
where we multiply each number by its location, the sum is now i x e, since 
the error occurred in position i. Since division is possible with these 
symbols, we can divide the second equation by the first, to give the 
position i on the right-hand side. But once we know the position i, we can 
correct the error. 


For example, 11... 1 is a word in the code. If we receive an 8 in position 4, 
instead of 1, then the size e of the error is 7, and the position of the error i 
is 4. From the message 1111811111111, we obtain the equations 


fo tt tito t+...t+%2=7, 
and 
11 + 2% +3r3t+...12r2 =2. 


So e = 7 andixe =2, and thus i = 2 + 7 = 28 + 7 = 4. Thus, an error of 7 
occurred in position 4, as expected. 
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The codes employed on compact discs work in a similar way, but they use 
binary numbers. For 3-bit numbers, we have eight possibilities which are 
added and multiplied using the following tables. Extending this idea to 8- 
bit numbers gives us the codes used for compact discs. 


If the disc is badly scratched, so that the errors occur in ‘bursts’, the 
compact disc system uses two codes, rather than one. The first encoder 
takes groups of 24 information symbols and produces groups of 28 
intermediate symbols. These are then fed through delay lines of 
increasing lengths. Then the second encoder converts groups of 28 symbols 
into groups of 32 output symbols. This means that two successive 
information symbols are now at least 32 symbols apart when encoded: the 
information is interleaved in time. These groups of 32 symbols are sent to 
the modulation encoder and then written to the disc. 


The decoding process is the reverse of the encoding process. Information 
read from the disc is passed through the modulation decoder, and then 
collected in groups of 32 code symbols. Each group is processed by a 
decoder, to give a group of 28 intermediate symbols. This decoder can 
correct up to two errors. If it detects more than two errors, it replaces the 
entire output by a group of special symbols to indicate a problem. The 
groups of intermediate symbols are then passed through delay lines of 
decreasing lengths, and another decoder, which can also correct up to two 
errors, processes them to give groups of 24 information symbols. Finally, 
these groups of 24 information symbols are turned into the analogue output 
signal. 


When there is a ‘burst error’, it affects several adjacent code symbols and 
the first decoder cannot correct them all. But it will almost certainly 
detect them, and its output will contain the special symbols. Since the 
delay lines spread these out, the second decoder only has to deal with one 
error in a group, and it can easily correct that. So, even if your copy of the 
disc is somewhat damaged, the result is a remarkably faithful 
reproduction of the original signal. 
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TV7 Transporting flowers 


Presenters: Roy Nelson, Andrew May 


Twenty-eight miles south-west of Land’s End lie the Isles of Scilly. 
Although St. Mary’s is the largest of the islands, it is only two-and-a- 
half miles long. Bathed in the warm waters of the Gulf Stream, the 
islands enjoy an equable climate which favours the growth of many sub- 
tropical flowers not found elsewhere in the British Isles. It is this which 
gives the islands one of their most important industries: the cultivation 
and sale of flowers. The growing season extends from November to May, 
and the farmers on Scilly frequently have flowers available for sale long 
before any of the mainland growers. Because the flowers are often picked 
in bud, the final customers get the maximum advantage from the flowering 
blooms. 


The centre of St. Mary’s is largely devoted to the growth of flowers from 
bulbs brought originally from the Mediterranean by the island’s seamen. 
Andrew May owns the Seaways Flower Farm, following in the footsteps of 
his predecessors who have been growing flowers in Scilly for over a 
hundred years. They grow varieties of narcissus, such as Avalanche, one of 
the older varieties, and the variety on which the industry was really 
founded, Grand Soleil D’Or. This latter variety accounts for about half of 
the acreage grown on the islands. 


The fields are small, protected from the wind by tall hedges of 
Pittasporum. However, these granite islands are exposed to strong gales 
from the Atlantic, and there is often a rush to gather in the flowers from 
the fields when they are threatened by storms at the crucial period when 
they are due to be picked. 


After picking, flowers are bunched in tens and packed fifty bunches to a 
box. The boxes are then stored in a cold-room which keeps them at a 
temperature of 2° Celsius. At this temperature the blooms do not 
deteriorate in less than six weeks. As many as 2000 boxes of flowers can be 
handled in one night — that’s one million blooms! 


To market their flowers, growers on the island work together in a 
cooperative. Orders are taken from shops and supermarkets, and each 
order gives rise to the problem: how should those flowers be transported to 
the customer? Once the flowers have been picked on the individual farms 
and packed into boxes, the growers have to distribute those flowers to 
their customers throughout the UK and Europe. 


The most difficult element of that is to get the boxes to Penzance over the 
twenty-eight miles of water. Once the boxes are in Penzance, the growers 
can solve the next stage of the problem quite easily. They get the flowers 
to Penzance either by sea in a freighter that takes about four hours to 
make the crossing, or by fixed wing aircraft, which takes about fifteen 
minutes to Land’s End, just west of Penzance, or by the helicopter which 
goes to Penzance itself. 


Although the sea freighter has the greatest capacity and can take the 
entire consignment with ease, the speed of air transport may be needed for 
urgent deliveries or important customers. The helicopter service is very 
reliable, and so is frequently used for part of the consignment. But the 
decision is not always straightforward. When the flowers arrive at the 
heliport, problems can often begin. The helicopter can carry up to 300 
boxes, but it may not be able to do so in certain weather conditions. If the 
atmospheric pressure is very low, then the lift of the helicopter can be 
reduced by as much as half a ton, and the growers may be told that not all 
of their flowers can be carried to the mainland. Since the largest 
consignment does not necessarily produce the greatest profits, an 
interesting mathematical problem arises. 
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The knapsack problem 


Let us suppose that there are four consignments that could be sent by 
helicopter: A, B, C and D, with the following weights and values. 


consignment A B 2 
weight (kg) GisriI; riivd 9 
value (£) a. & 7 @ 


On today’s flight, the helicopter can carry a maximum of 20 kg. It cannot 
manage all four consignments, so which ones should we select so as to 
maximize the total value taken? 

If we choose consignment B, with value £68, we cannot then include 
consignment D, because this would take us over the weight limit. A better 
choice would be to select consignments B and A, with a total value of £101. 
But can we do better? Let’s check all the cases. We'll represent each 
choice by a four-digit binary number; for example, choosing consignments A 
and C corresponds to the binary number (1, 0, 1, 0). 

The solution we found before was consignments A and B that is, (1, 1, 0, 0) 
with total value £101. Is this the best possible solution? 

To find out, we use a systematic method, adding items one at a time, and 
using a tree to show the ones we have included. 


We start with no items at all. 
(0, 0, 0, O)@ 


There are four choices for adding on an item. 
(1, 0, 0, 0) 
(0, 1, 0, 0) 
(0, 0, 1, 0) 
(0, 0, 0, 1) 


(0, 0, 0, 0) 


From (1, 0, 0, 0), we can add on any one of the following three choices: 


(1, 1, 0, 0) 
(1, 0, 0, 0) < (1, 0, 1, 0) 
(0, 1, 0, 0) (1, 0, 0, 1) 


(0, 0, 1, 0) 
(0, 0, 0, 1) 


(0, 0, 0, 0) 


But from (0, 1, 0, 0), only the choices (0, 1, 1, 0) and (0, 1, 0, 1) are allowed; 
we never add a 1 to the left of an existing 1. 

Note that, we do not need to check all sixteen possibilities. To see this, we 
check the weights and values. With no items, the weight and value are 
both zero. If we add the option (0, 1, 0, 0) with the greatest value (£68), 
we can branch in two possible ways, one of which exceeds the weight limit 
and so can be ignored. 


(1, 0, 0, 0) 

(0, 1, 0,0) — (0,1,1,0) weight 16, value 89 
(0, 0, 0, 0) anes (0,1,0,1) weight 21, value 116 

(0, 0, 1, 0) 

(0, 0, 0, 1) 


Similarly, (1, 1, 1, 0) has weight 22kg, which is also over the weight 
limit. So although we could add another item, that will also be over the 
limit, and so we do not have to consider it. 


This gives only a small saving in time, but with larger examples, the 
saving can be very significant indeed. You may like to check that the 
maximum value arises from the binary number (1, 0, 1, 1) — that is, items A, 
C and D with total weight 19 and total value £102. 
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So choices have to be made when we send goods on a route with limited 
capacity. And if we want to maximize the value of the goods carried, we 
can use the above method to work this out, without trying out all the 
possibilities. Details of this method were discussed in Graphs 4. 


Of course, the flowers that are not taken by air still have to be 
transported, and the alternative is sea transport. Although slower, this is 


less expensive, and there is no capacity limit. But in either case, air or | 


sea, this is just the first stage of the journey. 


All of the flowers exported from Scilly arrive in Penzance, and it is then 
of the utmost importance to get them to the customers without delay. Now 
the cost of transporting flowers to different customers varies according to 
two distinct factors: the mode of transport used, and the location of the 
customer. 


The transportation problem 


Thus a very interesting problem arises. Suppose that there are two orders: 
one from Birmingham, for 1100 boxes of flowers, and one from Edinburgh, 
for 600 boxes. There are three means of transport, each with a capacity 
limit: 200 boxes, by plane; 500 boxes, by lorry; and 1000 boxes, by train. The 
following table shows the cost per box for each mode of transport. 


plane _lorry train 


Birmingham 


Edinburgh 15 10 3 
To find the cheapest way of satisfying the customers’ requirements, we 
first find a feasible solution, ignoring the cost. We begin with the 
cheapest option, sending 1000 boxes to Birmingham by train. The other 100 
boxes at Birmingham are provided by lorry, as that is the next cheapest. 
The requirement at Birmingham is now satisfied, so we send 400 boxes to 
Edinburgh by lorry. Finally, we send the remaining 200 boxes to Edinburgh 
by plane. We can represent this situation by the following bipartite 
graph. 

(200) plane 


© Birmingham (1100) 


0 Edinburgh (600) 
(1000) train 


So that is a feasible solution, but is it optimal? When we calculate the 
total cost, we get £9500. But can we do better? Suppose we send one box by 
train to Edinburgh, instead of to Birmingham. We could then send one box 
by lorry to Birmingham, instead of to Edinburgh. If we calculate the total 
cost, we now get £9496 — a saving of £4. To gain the maximum saving in 
this way, we divert the remaining 399 boxes as well, and now the total 
cost is £7900 — a considerable saving. 


(200) plane 


© Birmingham (1100) 


0 Edinburgh (600) 
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The way we achieved that improvement was by altering the flow in the 
cycle on the left, interchanging lorry and train. Now let us try the cycle on 
the right, interchanging plane and train. 


plane 
Birmingham Birmingham 
lorry 
Edinburgh : Edinburgh 


train train 


Diverting one box saves £3 on the plane, and costs an extra £1 on the train, 
with a net saving of £2. So we can divert 200 boxes, giving a total cost of 
only £7500 — again a useful saving. 


(200) plane 


0 Birmingham (1100) 


0 Edinburgh (600) 
(1000) train 


The only cycle that we have not considered interchanges plane and lorry. 
But now, diverting one box saves £3 on the plane, but costs an extra £5 on 
the lorry, so this does not give us a further saving. We have, in fact, 
arrived at the optimal solution. 


So that is another problem that can be solved efficiently. Where there 
are different methods of transport available, with different costs and 
capacities for different destinations, we can again reach the best solution 
without trying out all the possibilities. An algorithm for solving any 
problem of this type was given in Networks 3. 


But there are other considerations, too. There is a trade-off, because the 
further you distribute your narcissi, which are a relatively low-value 
item in the cut-flower world (when compared to orchids, for example), 
then the higher the proportion the transport cost is of the overall value. 
This has to be taken into consideration, and most of the flowers in fact go 
to the Midlands and to London to be processed and sold to supermarkets 
and florists. However, the Scillonian growers also supply Scotland on a 
regular basis, and send to Holland for onward distribution throughout 
Europe. 


Holland is a good place to send flowers for auction, because nearly 60% of 
the entire world export market is based there. The largest of these 
markets is at Aalsmeer. Flower growers from all over the world send their 
flowers to the auction rooms there. 


On a typical day there can be nearly 40 million blooms in the vast hall at 
Aalsmeer. With a floor area of over 700,000 square metres, it is one of the 
largest buildings in the world, and as well as the storage areas, it houses 
five separate auction rooms. Every day, around 15 million flowers and 2 
million plants are sold there. There is a great advantage in bringing such 
a large share of the world’s flower supply into a single place. Market 
forces ensure that all the growers get fair prices for their products. 


Before sale, the flowers are loaded onto trolleys — an average of ten 
thousand a day. Each trolley is then guided along a chain track into one of 
the auction rooms and pauses below one of the clocks. Inside the auction 
room, selling takes place by the traditional ‘Dutch Auction’ method: the 
auctioneer starts with a high price and works downwards until somebody 
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bids. But modern technology is used, too. Instead of the price being shouted 
out, it is shown on computer-controlled clocks with lights that start at 100 
and move rapidly around the dial down to 1. The unit may be anything 
from a cent to a guilder. 


Each buyer presses a button when the clock gets down to a price he or she is 
willing to pay for the current batch on offer. Each registered buyer has an 
address in the huge building, and everything bought by that buyer jis 
sorted out and delivered on a trolley to that address. The ‘distributors’ are 
essential to the whole operation of the auction; they link the trolleys into 
long trains and shuttle them all over the building. 


Once the flowers have been paid for, they can be loaded into vans and 
lorries for despatch. Some buyers come from small concerns, maybe with 
just a single shop in a nearby town, and might need just one or two small 
lorries. But there are also buyers from large international supermarket 
chains. As they want to get the flowers to Schipol airport as quickly as 
possible, they use several of the largest trucks available. But there are 
still choices to be made. 


The bin-packing problem 


Having bought all the flowers needed, the buyers may sometimes find 
themselves faced with a problem. If there are too many flowers to fit into 
one lorry, how many lorries are needed? This is an example of the bin- 
packing problem. Suppose that we have a collection of ten items of 
different sizes, and several bins (or lorries) each with a capacity of 20 
units, how should we pack the items so as to use the smallest number of 
lorries? 


item: A 2 ©€.2. €. © -o oe wes 
size: 4 9 9 4 7 5 6 = B11 


A basic approach (the first-fit packing algorithm) takes the items in the 
order they occur, and packs each item in the first lorry with enough spare 
capacity to take it. (As we saw in Networks 2, such a method is called an 
on-line algorithm.) So the first item goes in lorry 1. The second item will 
also fit in lorry 1; the third item will not fit in lorry 1, so we put it in lorry 
2;andsoon.... 


20 


lorries 


Altogether, that uses six lorries. But we can do better. It makes sense to 
pack the largest items first, and then to fit the smaller ones into the gaps. 
(Such a method, where we can rearrange the items, is called an off-line 
algorithm.) So we start by rearranging the items in decreasing order of size. 


item: eee eee es * 
size: Se ee ee, ee ee eo ee oo ee 
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The first three big items go in lorries 1, 2 and 3, and then the gaps are 
filled up. 


20 


lorries 


So this method (the first-fit decreasing packing algorithm) uses only five 
lorries, which is fewer than before. But notice that the final lorry is 
nearly empty. Is there a way of juggling the items to get them into just four 
lorries? In fact there is — like this. 


20 


lorries 


But we might have to discover that by trial and error. For the bin-packing 
problem, there is no known method that is always guaranteed to give the 
best answer. 


So we have seen three different mathematical problems associated with 
flowers, and we have suggested methods for solving them. For the 
knapsack problem of loading the helicopter, we can always find the 
optimum solution, although there is no known polynomial-time 
algorithm. For the transportation problem, we can find the best solution in 
polynomial time. But the bin-packing problem is different. Although we 
have suggested two simple methods, they rarely give the optimum 
solution in practice. This has turned out to be a problem that 
mathematicians have not yet solved, except by an exhaustive method. 


So even a simple bunch of flowers can involve much of mathematical 
interest. 
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